The total nucleon-nucleon cross section at large N c 
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It is shown that at sufficiently large N c for incident momenta which are much larger than the 
QCD, the total nucleon-nucleon cross section is independent of incident momentum and given by 
^totai _ 27rlog 2 (7V c )/(m4) . This result is valid in the extreme large N c regime of log(iV c ) > 1 and 
has corrections of relative order log (log(iV c )) / log(iV c ). A possible connection of this result to the 
Froissart-Martin bound is discussed. 
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The large N c limit of QCD and the 1/N C expansion 
have been of great interest since introduced by 'tHooft 
nearly 40 years ago [l[. While the approach to date has 
not provided a path by which quantities can be calcu- 
lated ah initio directly from QCD except in special cases 
such as QCD in 1+1 dimensions |2[ or QCD in the limit 
of heavy quark masses [3j, |j| , it has provided a qualitative 
understanding of many aspects of hadronic phenomena. 
Witten's extension of the analysis to include baryons has 
played a critical role [3] . One of the remarkable features 
of baryons is the emergence of a contract SU(2Nf) sym- 
metry at large N c [5( which has allowed predictions of 
both ground state baryons [5|] and excited baryonic res- 
onance j6| . In this work we will focus on the nonstrange 
sector and assume exact isospin invariance. 



Implications of large N c QCD for nuclear physics were 
first explored in Witten's seminal paper on large N c 
baryons Q. A key result of this analysis is that the 
nucleon-nucleon interaction has a strength which scales 
as Nl and a range which scales as N®. Moreover, 
nucleon-nucleon scattering with fixed incident momen- 
tum has no smooth large N c limit. However, a sensible 
time-dependent mean-field description emerges if the ini- 
tial velocity is held fixed at large N c (that is, that mo- 
mentum scales linearly with N c given that the mass is 
linear in N c ). While there has been significant work on 
various aspects of nuclear physics at large N c , such as 
treatments of the NN potential [7|, the phenomenologi- 
cal relevance of the large N c limit for nuclear physics is 
far less clear than for hadronic physics Q. Despite this, 
it is of interest to understand the N c scaling behavior of 
quantities of interest in nuclear physics. One quantity 
which has received comparatively little attention except 
for a recent paper on its spin flavor dependence [9| is the 
total nucleon-nucleon cross section (with the effects of 
electromagnetic interactions removed). This is unfortu- 
nate since, as will be shown in this letter, the total cross 
section is truly remarkable in that it can be computed 



analytically when N c is sufficiently large: 



total 



27rlog 2 (iV c ) 



»K 



(1) 



Equation (Q]) holds for all spin-isopsin channels in the 
regime where the incident momentum is much larger 
than Aqcd! corrections to Eq. ([T]) are of relative or- 
der log (\og(N c )) / \og(N c ) . Formally N c needs to be ex- 
tremely large for Eq. (JT]) to hold and it is not obvious that 
the result is phenomenologically relevant for the physi- 
cal world of N c — 3. In any event, the result is of real 
interest from the perspective of theory. 

To gain insight, it is useful to first consider a simplified 
"toy problem" of scattering of two nonrelativistic spin- 
less particles of mass M interacting via a central poten- 
tial which falls off exponentially at large distances. The 
problem has a control parameter, A, which controls the 
scaling of the potential strength (but not its range) , the 
mass, and the initial relative momentum: 



M = XM V(r) = XV (r) k = Xk: 



(2) 



where the quantities with a tilde are independent of A. 
The classical trajectory followed by a particle in this 
problem depends on the impact parameter b and p, V(r), 
and M but not on A: if the quantum scattering is de- 
scribed well classically, then the differential cross section 
will be independent of A. The parameter, A, however, 
controls the region of validity of a semi-classical descrip- 
tion in an underlying quantum scattering problem; the 
classical limit corresponds to large A. Of course, by de- 
sign this problem mirrors the N c scaling rules of NN scat- 
tering with A playing the role of N c . 

The potential is central and the partial waves are in- 
dependent. Using the scaling rules in Eq. © and the 
Schrodinger equation for a given partial wave, the phase 
shifts can be shown to scale as 



5 l (k) = 6 xl (Xk) = XS t (k), 



(3) 



'Electronic address: cohcn@physics.umd.edu 



where / = l/X is introduced for convenience and S is in- 
dependent of A. Corrections are of relative order 1/A. To 
see this, parameterize the radial wave function for a given 
partial wave as the product of a phase and an amplitude, 



"0i( r ) = e 1 *'^^^)], an( i take I to be proportional to 
A. It is easy to show self-consistently that at large A, <& 
is proportional to A while |V'i( r )l is slowly varying and 
independent of A. This is precisely what one expects if A 
acts as the control parameter for the semi-classical limit. 
The total cross section for central potentials is given by 



total 
J toy 



(k) 



4tt 



^(2; + l)sin 2 (^(fc)) 



4tt 



(4) 



— / dl 21 sin 2 (XS^k)) 



where Si is the phase shift for the I th partial wave and 
the integral expression becomes exact as A^oo. 

In the integral of Eq. (QJ, focus on the range from It 
to li . For large A, sin 2 oscillates rapidly and averages to 
\ (up to corrections of order 1/A) over this region; the 
contribution to the cross section becomes 



A< al (fc) = 2tt(^ b{) 



. , , / I 
with b — — = — . (5) 
k k 



shifts are well approximated [ll| by 



l/k 



k 2 - P/r 2 



zdr . 



(7) 



where jl is the reduced mass divided by A. 

For concreteness, take the form of the potential 
to be the sum of Yukawa interactions: V(r) = 

J2n C n c ~r wnere the C n are strength parameters 

independent of A and the r„ are the ranges. Note that at 
large A, 6 cut is large and the integral in Eq. ([7]) is domi- 
nated by the longest-range contribution to the potential. 
Evaluating the integral yields 



CM 
^~-K (l/(kr ) 



.^Koibi/ro)) . (8) 



where r$ is the longest range in the potential and Co 
is the associated strength. For large values of bj/ro it 
is legitimate to use the asymptotic form of the Bessel 
function when inverting this relation; doing this yields 



This is twice the geometric cross section associated with 
impact parameters from b\ and bi ; the factor of 2 is due to 
a nearly forward diffractive scattering contribution equal 
to the geometrical contribution [10]. At infinite A there 
is no bound on the Is that contribute implying that the 
total cross section diverges as A — > oo. 

The quantum cross section is finite because the phase 
shifts approach zero as I — > oo. For any finite value 
of A, there is a regime of sufficiently large I such that 
8 ~ A -1 and the sin 2 term does not oscillate rapidly to 
yield an average of i. The phase shift in this regime 
rapidly becomes small and makes small contributions to 
the total cross section. The value of I beyond which the 
rapid oscillations effectively turns off depends on A and 
gets pushed off to infinity as A —¥ oo. 

To see this, start with the integral form of Eq. ((4]) and 
change variables into an integral over 5. It is a simple 
matter to show that 



total 
T toy 



2irbi ut +2nk- 



dsl^L) sin 2 (A5)+0(A- 1 ) 



(6) 

,72 - 

where ^- is treated as a function of S in the integral and 
6 cn t is an arbitrary but fixed "cutoff" phase shift of order 
A ; b cu t is the impact parameter associated with <5 cu t. 
The arbitrariness in the choice of (5 cu t is compensated 
by the integral in Eq. ([6]). If the integral in Eq. ([6]) is 
parametrically smaller than 2nb 2 ut when (5 cu t is of order 
A , then up to parametrically small corrections <r'°* al = 

As shall be shown self-consistently, the total cross sec- 
tion is dominated by the behavior for large impact pa- 
rameters, or equivalently, large I. In this regime, each 
partial wave is semi-classical and the potential is much 
smaller than the centrifugal barrier. Thus, the phase 



' 2 ^ Sf# , 



(9) 



where W is the Lambert function. As x gets very 
large W(x) — > log(x) (reflecting the dominantly expo- 
nential behavior of Kq) with corrections of relative order 
log(log(x)) /log(x). At large I, the phase shifts become 
small; the log is dominated by 6t; bi ~ — rolog(Si) and 
frcut = — ?"o log(<5 cu t/A). Thus up to corrections of order 
A 



bent = r log(A) 



(10) 



and at large A is o-to y al = 27rr 2 log 2 (A) provided the in- 
tegral in Eq. ([6]) is parametrically small — which, as will 
be shown shortly, it is. 

Note that the sensitivity to the particle's mass and to 
the strength of the potential are contained in the order 
A correction terms to Eq. (ITU1) . Note further, that the 
sensitivity to the choice of S cu t is also contained in the 
A correction terms. Since the dependence of the choice 
of 8 cu t is compensated by the integral in Eq. ©, it fol- 
lows that the integral is also parametrically of order A 
and makes a negligible contribution to the cross section 
at large A. Thus, al°y al — 27rr 2 log 2 (A) with corrections 
of relative order log (log (x)) /log(x). With the substi- 
tutions tq — > l/m^ and A — > N c , this is of the form of 
Eq. (fTJ). It should be apparent that any power law pref- 
actor to the Yukawa potentials cannot alter this result at 
leading order. 

The result also holds for a relativistic version of the 
toy problem. Consider potential scattering in a relativis- 
tic two-body model (which lacks micro causality but can 
be consistently formulated as a quantum theory Il2l|). 
The basic set up remains intact: the partial wave decom- 
position still holds and Eqs. (Ql and (JU) remain valid. In 



the semi-classical regime with sufficiently large I, one can 
always cast the phase shift into the form of Eq. (JTJ) with 
V(r) replaced by (an energy dependent) V c e(r) whose 
form depends on the transformation properties of the in- 
teraction. If the longest range interaction in the model 
transforms as a Lorentz scalar (as in QCD), then at long 
range V c s(r) — V s (r)E/p,. Relativity affects & C ut only by 
renormalizing the strength of the longest-range interac- 
tion by an energy-dependent amount independent of A. 
Since the leading behavior of the total cross section does 
not depend on the strength, the relativistic toy model 
also has: 0"*°*^ = 2nr 2 log 2 (A) which corresponds to 
Eq. (HJ. 

Nucleon-nucleon scattering in large N c QCD is clearly 
more complicated than in the toy problem for several 
reasons: i) The nucleons have spin and the partial wave 
expansion for elastic scattering is necessarily of a coupled 
channel form; ii) there is an emergent spin-isospin sym- 
metry at large N c i5J; as a result of this symmetry the 
A and a whole tower of baryons are stable and nearly 
degenerate with the nucleon at large N c . The emergent 
symmetry implies correlations between channels in scat- 
tering p], [3, LLJ] ; hi) There are inelastic channels due to 
meson production. However, as discussed below, the re- 
sult in Eq. ((T|) is quite robust and is unaltered by these 
complications. 

To treat the full problem, the generically strong (or- 
der N c ) nature of the NN interaction must be encoded 
in a model-independent way. The potential model treat- 
ments of the toy problem are inappropriate to this prob- 
lem and, in any event, the potential is intrinsically un- 
physical which can lead to subtleties in the N c counting 
[13]. The physically relevant object is the S-matrix for 
nucleon-nucleon scattering, S NN . Its elements for elastic 
scattering can be denoted Sf^.y a , where a (a') specifies 
the spin and isospin configuration of the incident (final) 
states of the two nucleons. Conservation of angular mo- 
mentum and isospin along with the fermion nature of 
nucleons constrain the form of S NN ; for example, the 
matrix elements are zero unless I' = I — 1,1,1 + 1. 

The strong NN interaction at large N c cannot be en- 
coded by having S NN scale linearly in N c ; it is bounded 
due to the unitarity. To proceed, note that in the toy 
problem, Eq. ^ the phase shift, i.e the logarithm of the 
S matrix in the partial wave channel, scales with JV C . As 
will be discussed in detail in a future publication, this 
behavior should hold generically for diagonal matrix ele- 
ments of the S matrix in large N c QCD. 



log(S™ )^2^, a 



2<a 



Ka 



N r 



(11) 



where the S-matrix is for nucleon-nucleon elastic scatter- 
ing at fixed initial velocity. Note that 8i^ a is not real in 
general; the imaginary part reflects scattering out of the 
original channel either to other elastic channels (with dif- 
ferent final I or a) or to inelastic channels. Both the real 
and imaginary parts are expected to scale with N c . This 
applies to all physical channels (eg., two neutrons with 
spins aligned with the beam). 



The total cross section in multi-channel problems with 
the initial nucleons in spin-isospin configuration a can be 
expressed in terms of the diagonal matrix elements of the 
S matrix: 



total 



(k) 



(12) 



% J2(2l + 1) (1 - cxp(-2^ a (fc)) cos(2«J* (fc))) 



k 2 ^ 
1 



4-7T 
I 2 



d? (l - exp (-2N c 6f a {N c k)\ cos (2N c Sf a {N c k) 



where k = N c k ; the first form is general [ll| and the 
second form builds in the N c scaling of the phase shifts. 
The integral form becomes exact in the limit N c — >• 00. 

Note that Eq. {T5J coincides with Eq. Q of the toy 
problem if one sets 8j a to zero. Moreover, Eq. ([5]) con- 
tinues to hold providing the integrand is in the regime 
where either the real or imaginary parts of 8 (or both) 
are of order unity; in the case of the real part it is due 
to rapid oscillations as in the toy problem; in the case of 
the imaginary part it holds due to an exponential sup- 
pression. As in the toy problem, the total cross section 
is determined by where 8 ceases to be of order unity and 
becomes of order 1/N C . It is easy to see that if, as a 
function of I, 8{ a approaches zero at least as rapidly as 
8^ a , then the total cross section at leading order will be 
determined by where 8- 1 drops to order of 1/N C . 

It is clear that 8{ a does approaches zero more rapidly 

than 8f- a : consider going to sufficiently large I so that the 
phase shifts are accurately described by the Born approx- 
imation for the longest range part of the interaction — 
one-pion exchange. In that regime the 8^ a is small but 
non-vanishing while 6j a vanishes at the first Born ap- 
proximation level only arising at second order. Moreover, 
it is clear that at very large I the real part of the phase 
shift is dominated by the Born approximation contribu- 
tion to one-pion exchange which drops off exponentially 
in exactly the same way that it does in the toy model. 
Accordingly the result in the toy model carries across and 
Eq. |T]) follows exactly as in the toy model. Note that this 
exponential fall-off holds for any physical spin and isospin 
configuration of the initial baryons; pion exchange domi- 
nates regardless of the initial spin orientations or whether 
the two nucleons are the same or different. Thus Eq. ([T]) 
also holds for any initial configuration and the leading 
order cross section is spin and isospin independent. This 
is consistent with the analysis of Ref. ]9], although it 
is more restrictive than the most general leading-order 
result deduced there. 

It is worth observing that at N c , elastic scattering 
should account for 1/2 of the total scattering or more. 
In the regime of I where Sj a ~ N c , the contribution to 
the scattering looks like a black disk for which diffractive 
scattering is 50%. There may, in principle, also be sub- 
stantial contributions for the regime where 8{ a is small 
but 81 a ~ N c ; such contributions will have elastic contri- 
butions of greater than 50%. 



The log 2 form of the cross section in Eq.fi} is strik- 
ingly similar to the Froissart-Martin bound [l5| . At large 
Mandelstam s, considerations of unitarity, analyticity 
plus the knowledge that the pion is the lightest excita- 
tion in the system serve to bound the growth of the total 
cross section at large s 



total ^ 



11 1 2 

— lo s , 
mi \ s 



(13) 



where sq is a reference scale. This similarity may not 
be accidental. Note that the natural regime for nucleon- 
nucleon scattering at large N c is for fixed velocity [3] 
which in turn implies that s ~ A 2 . Provided that s 
does not also scale with N c , the bound becomes a totaX < 

-tt log 2 ( ^j^- ) where s is independent of N c . If one takes 

the large N c limit prior to the large s limit, and keeps only 
the leading behavior, one has a total < ■£&■ log 2 (N c ) up to 
corrections of relative order 1/ log(iV c ). Note that Eq. ^ 
satisfies this inequality by exactly a factor of i. This 
factor of \ is suggestive. The derivation of the Froissart- 
Martin bound requires unitarity. However, if one looks 
at the integral form of Eq. (|12[) it is clear that in the 
region of dominant contribution, the integrand at large 
N c is precisely h of its unitarity bound (which occurs at 
S R = 7r/2, S 1 = 0). Thus, the present result is natural in 
light of the Froissart-Martin bound. 

To what extent is this result applicable to the physical 
world of iV c =3? In the physical world, the total cross 
section for a/s well above AQCD is approximately inde- 
pendent of s[18} as would be expected from the large N c 
analysis. Over three orders of magnitude in ^fs, from 
1.5GeV < yfs < 1200GeV, the cross section for proton- 
proton scattering varies by only about 25%. Moreover, 



the cross section is dominantly spin and isospin indepen- 
dent [1, [l6| as predicted to occur at large 7V C . These 
results may suggest that the large N c analysis is of phe- 
nomenological relevance for the physical world of N c = 3. 
However, this is not clear. For example, above 1.5 GeV, 
the total cross section is predominantly inelastic; the elas- 
tic cross section is typically less than 4 of the total cross 
section and by t/s of several 10s of GeV, it drops to under 
20%. While it has been argued that at truly asymptoti- 
cally high energies [17( it approaches \, at large N c this 
behavior is expected for all s well above J^-\ rm QCD- This 
implies that a substantial part of the cross section comes 
from regions where both the real and imaginary parts 
of the phase shift are small. This is at odds with the 
behavior expected at large N c where, as noted above, 
elastic scattering should be 50% or higher. Given this, 
it is perhaps not too surprising that the absolute predic- 
tion of Eq. ((TJ) that cr total « 150mb is significantly larger 
than the empirical value of approximately 40 mb. Ul- 
timately, the reason that the large N c analysis for the 
magnitude of the total cross sections is not very predic- 
tive for the N c = 3 world is quite understandable. For- 
mally, one expects the analysis to be predictive only when 
log(N c ) > 1. Clearly this is not true for N c = 3. What- 
ever, the phenomenological significance for the world of 
N c = 3, the fact that at large N c the total cross section 
is calculable is, at the very least, of theoretical interest. 
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